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1.  Summary 

Processes  related  to  the  supply,  mixing,  and  removal  of  ocean  bubbles  and  their 
impact  on  low-frequency  (0-1000  Hz)  acoustic  surface  scattering  strength  (SSS)  in  the 
upper  ocean  have  been  investigated.  A  key  finding  is  that  gross  site-to-site  differences 
between  Critical  Sea  Test  (CST)  and  Chapman-Harris  SSS  observations  are  explained  by 
the  mean  ocean  temperature,  biological  productivity,  and  wave  conditions  at  each  site. 
Furthermore,  wave  energy  dissipation  rate,  estimated  from  wave  spectra  using  wave 
energy  flux  considerations,  was  found  to  be  a  better  descriptor  for  ocean  whitecaps, 
related  to  upper  ocean  bubble  supply,  than  wind  speed.  These  results  demonstrate  that  a 
model  could  be  developed  for  estimating  SSS  statistics  in  coastal  areas  using  readily 
available  (by  satellite  and  operational  models)  environmental  inputs. 

2.  Background 

As  an  important  element  of  imderwater  sonar  performance,  SSS  measurements 
have  been  made  in  a  variety  of  oceanographic  conditions  by  numerous  investigators. 
The  recent  series  of  CST  experiments  (1988-1992)  has  provided  a  rich  set  of  observations, 
obtained  with  a  consistent  technique,  in  six  different  environments.  The  calibrated  CST 
results  can  be  readily  compared  with  earlier  measurements  such  as  those  of  Chapman 
and  Harris  obtained  North  of  Bermuda  in  1962  {Reference  1). 

Surface  scattering  results  from  the  CST  program  have  been  reviewed  in  a  series  of 
reports  by  M.  Nicholas,  P.  Ogden,  F.  Erskine,  and  R.  Gauss  of  NRL  and  F.  Henyey  and 
E.  Thorsos  of  APL/UW.  The  CST  results  most  relevant  to  our  study  are  listed  here: 

•  SSS  is  strongly  linked  to  short-term  (~  1  h)  wind  history. 

•  Supporting  environmental  data  suggest  that  SSS  is  more  closely  associated 
with  tenuous  bubble  clouds  than  with  breaking  wave  events. 

•  Acoustic  model  comparisons  with  CST  results  indicate  that  near-surface 
(~  1-3  m)  microbubbles,  entrained  by  wave  orbital  motions,  contribute 
most  to  SSS. 

•  Significant  differences  (6-7  dB)  are  noted  between  CST-4  and  CST-7  SSS 
results  at  similar  wind-forcing  conditions  (see  Figure  1).  Unexplained  dif¬ 
ferences  also  exist  between  other  CST  data  sets  and  with  Chapman-Harris. 

•  CST  1,  2,  3,  4,  5,  and  7  results  used  to  construct  NRL  Ogden-Nicholas- 
Erskine  (ONE)  empirical  SSS  model.  Environmental  input  is  the  1-hr 
backaveraged  wind  speed. 
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Actual  values  between  these  sites  differ  by  as  much  as  8  dB.  Adapted  from  Nicholas, 
M.,  P.  M.  Ogden,  and  F.  T.  Erskine,  Improved  Empirical  Descriptions  for  Acoustic 
Surface  Backscatter  in  the  Ocean,  submitted  to  lEEEJ.  Oce.  Engineering. 

Figure  1  Scattering  Strength  Differences  Between  CST-4  and  CST-7 


A  crucial  issue  remaining  from  the  CSX  program  is  explanation  of  the  large  site-to-site 
differences  in  SSS.  It  was  anticipated  that  results  from  the  Acoustic  Surface  Reverbera¬ 
tion  Experiment  (ASREX),  conducted  west  of  Bermuda  over  a  3-month  period  during 
the  1993-1994  winter  storm  season,  would  identify  the  environmental  factors  missing 
from  the  Chapman-Harris  and  ONE  SSS  models.  Unfortunately,  calibrated  SSS  levels 
from  ASREX  are  not  yet  available  for  direct  comparison  with  the  CSX  results.  How¬ 
ever,  the  extended  analysis  of  CSX  observations  reported  here  has  led  to  a  possible 
explanation  for  site-to-site  differences  in  SSS. 
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3.  Research  Objectives 

Our  h5^othesis  is  that  site-to-site  variations  in  SSS  are  a  result  of  physical  and  bio¬ 
logical  factors  related  to  the  supply,  mixing,  and  removal  of  bubbles  in  the  upper  ocean. 
As  the  wind  speed  input  of  SSS  models  appears  to  adequately  prescribe  the  local,  short¬ 
term  forcing  of  the  upper  ocean,  the  missing  environmental  factors  must  relate  to  the 
background  setting  of  each  site  (swell,  mixing,  temperatures,  dissolved  gasses,  surfac¬ 
tants,  etc.).  Hence,  our  primary  research  objectives  were  the  following: 

(1)  Examine  the  influence  of  bubble-related  environmental  factors  on  SSS  with 
particular  emphasis  on  parameters  linked  to  seasonal  and  geographic  vari¬ 
ations  in  near-surface  bubble  populations.  The  need  to  identify  the  'miss¬ 
ing  link'  for  explaining  site-to-site  differences  in  SSS  is  recognized. 

(2)  Determine  if  surface  wave  descriptors  can  be  used  to  augment  or  replace 
wind  speed  in  air-sea  process  models  such  as  those  for  air-sea  gas  flux  (by 
bubble  injection)  and  acoustic  surface  scatter. 

4.  Site-to-Site  SSS  Variability:  An  Explanation 
4.1  SSS  Prediction  Errors 

The  NRL  ONE  model  predicts  SSS  as  a  multiparameter  function  of  frequency, 
grazing  angle,  and  wind  speed.  Observation  departures  from  ONE  predictions,  aside 
from  measurement  errors,  result  from  environmental  influences  not  well  represented  by 
wind  speed  alone.  Here  we  isolate  these  departures  by  calculating  the  SSS  prediction 
error 


SSSe  -  SSSpi^g-SSS^i^^  , 

with  the  'ONE'  and  'obs'  subscripts  referring  to  model  predictions  and  observations, 
respectively. 

The  mean  SSSe  at  500  Hz  for  each  of  the  CST  and  Chapman-Harris  sites  appears 
in  Figure  2.  The  mean  levels  were  calculated  by  averaging  the  prediction  errors  for  each 
experiment  over  all  wind  speeds,  all  grazing  angles  >  10°,  and  over  the  frequency  band 
400-600  Hz.  Error  bars  depict  one  standard  deviation  about  the  means.  It  is  observed 
that  CST-2  and  CST-5  are  over-predicted  and  that  CST-1,  CST-3,  CST-4,  and  Chapman- 
Harris  are  under-predicted.  The  CST-7  values  are  close  to  zero;  this  is  expected  as  CST-7 
has  contributed  by  far  the  most  observations  used  for  the  ONE  model  fits.  The  differ¬ 
ences  between  CST-4  and  CST-7  levels  have  been  particularly  frustrating  to  the  research 
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CST-1  CST-2  CST-3  CST-4  CST-5  CST-7  Chapman- 


Observation  Sites 


Harris 


One  standard  deviation  error  bars  are  depicted. 

Figure  2  Mean  Surface  Scatter  Strength  Prediction  Error  {SSSqi^e  ~  SSSgi,s)  at  500  Hz 

community  as  both  experiments  occurred  in  the  Gulf  of  Alaska  under  similar  wind  and 
wave  conditions  (see  Figure  1).  Note  that  much  larger  differences  exist  between  the 
remaining  CST  experiments,  such  as  extreme  differences  between  CST-3  and  CST-5  of 
nearly  15  dB. 

Our  approach  is  to  empirically  relate  SSSe,  averaged  over  various  frequency,  graz¬ 
ing  angle,  and  wind  speed  regimes,  with  the  bubble-related  environmental  descriptors 
described  in  Section  4.2. 
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4.2  Environmental  Descriptors 

The  ONE  model  wind  speed  input  represents  those  environmental  processes  that 
influence  surface  scatter  on  short  (hourly)  time  scales  including  surface  wave  develop¬ 
ment  and  the  supply  of  bubbles  by  breaking  waves.  Missing  from  this  model,  however, 
are  environmental  descriptors  for  the  background  setting  of  each  site,  such  as  swell 
activity,  mixed-layer  temperature  and  depth,  gas  saturation,  and  surfactants.  Both  labo¬ 
ratory  experiments  and  theoretical  calculations  have  shown  the  importance  of  these 
background  conditions  on  processes  that  supply,  entrain,  and  remove  bubbles  from  the 
upper  ocean.  The  full  complement  of  supporting  environmental  data  collected  during 
CST,  combined  with  remote  sensing  products,  has  allowed  a  reasonable  test  of  their 
importance  to  SSS. 

The  environmental  descriptors  found  to  be  most  important  for  site-to-site  SSS 
variability  appear  in  Table  1.  Four  factors  are  employed  to  describe  bubble  supply, 
entrainment,  and  removal  processes  in  the  ocean:  wind  speed  {U^Q),  ocean  temperature 
(ssf),  significant  wave  height  {h^),  and  chlorophyll  concentration  (chi).  The  importance 
of  wind  speed  has  already  been  demonstrated  by  its  success  in  the  various  SSS  model 
formulations  (for  example,  Chapman-Harris  and  ONE). 


Table  1 

Environmental  Descriptors  for  Surface  Scatter  by  Bubbles 


Bubble  Issue 

Process 

Descriptor 

Supply 

Wave  breaking 

Wind  speed 

Ocean  temperature 

Entrainment 

Wave  mixing 

Wave  height* 

Removal 

Gas  dissolution 

Ocean  temperature 

Chlorophyll  concentration 

*  Mean  over  entire  test  duration. 


Significant  wave  height,  averaged  over  the  entire  test  duration,  is  employed  here 
to  represent  the  test  site  background  energy  setting  and  helps  account  for  the  effects  of 
swell  on  wave  breaking  and  mixing.  It  is  expected  that  high-energy  environments  will 
support  denser  and  deeper  bubble  populations. 

Sea  surface  temperature  (ssf)  significantly  influences  both  the  supply  and  removal 
of  bubbles  in  the  upper  ocean.  Both  surface  tension  and  viscosity  are  highest  in  cold 
water;  this  fact  offers  a  preliminary  explanation  for  reduced  whitecap  coverage  in  high 
latitudes  {Reference  2).  Water  temperature  also  controls  gas  solubility.  As  less  air  can 
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dissolve  in  warmer  water,  we  expect  longer  bubble  lifetimes  and  hence,  denser  bubble 
populations  with  water  temperature  increases.  CST  surface  temperatures  were  collected 
from  0.5-1.5  m  depth  and  represent  the  typical  ocean  temperature  within  the  near-sur- 
face  bubble  layer. 

Phytoplankton  biomass,  characterized  by  remotely  sensed  chlorophyll  concentra¬ 
tions,  will  also  influence  bubble  lifetimes  through  two  complementary  mechanisms. 
First,  surfactants  secreted  by  biological  communities  are  known  to  provide  a  stabilizing 
coating  to  seawater  bubbles.  Second,  phytoplankton  blooms  can  supersaturate  the  water 
with  respect  to  oxygen  and  other  atmospheric  gases.  Both  of  these  effects  will  decrease 
the  rate  of  bubble  dissolution  and  promote  larger  bubble  populations  in  biologically 
productive  areas.  This  effect  will  be  most  important  in  nearshore  and  coastal  waters 
where  biological  productivity  is  greatest. 

Monthly  mean  phytoplankton  pigment  (chlorophyll)  concentrations  for  the  CST 
and  Chapman-Harris  sites  were  obtained  from  the  Coastal  Zone  Color  Scanner  (CZCS) 
mission  results  on  the  World  Wide  Web  (http://seawifs.gsfc.nasa.gov/SEAWIFS.html).  This 
observation  set  covers  7.5  yr  from  October  1978  through  June  1986.  As  the  CST  program 
did  not  begin  imtil  1988,  the  CZCS  pigment  concentrations  can  only  be  used  to  represent 
tj^ical  conditions  at  each  site.  The  North  Atlantic  Ocean  images  for  the  CST-1,  CST-2, 
CST-3,  and  Chapman-Harris  sites  appear  in  Figure  3.  Black  areas  on  the  images  corre¬ 
spond  to  regions  with  no  data,  primarily  due  to  ice  and  cloud  cover.  A  dramatic  increase 
in  pigment  concentration  is  observed  at  high  latitudes  between  early  spring  and  late 
summer,  driven  primarily  by  the  availability  of  nutrients  at  high  latitudes  and  the 
increase  in  light  during  the  summer.  Note  in  particular  the  large  contrast  between  the 
CST-1  and  CST-3  sites,  and  the  highly  variable  nature  of  the  CST-2  region.  The  North 
Pacific  Ocean  images  for  February  (CST-7)  and  April  (CST-4)  appear  in  Figure  4.  The 

effect  of  spring  blooms  in  April  is  apparent  with  an  increased  level  of  productivity  at  the 
CST-4  site. 


4.3  A  Simple  Check  of  Hypothesis 

The  dependence  of  gross  site-to-site  differences  in  surface  scattering  on  the  bub¬ 
ble-related  environmental  descriptors  (Table  1)  can  be  demonstrated  with  a  few  simple 
empirical  tests.  SSSe  values  from  CST  and  the  Chapman-Harris  model  were  first  aver¬ 
aged  over  three  distinct  frequency  bands:  0-150  Hz,  400-600  Hz,  and  800-1000  Hz.  The 
data  from  all  runs  (individual  observation  sets)  within  each  CST  experiment  were 
included  in  the  averages  as  well  as  observations  at  all  grazing  angles  >  10°.  This  gross 
averaging  process  was  performed  to  minimize  uncertainty  due  to  random  measurement 
errors,  environmental  patchiness,  short-term  variability,  etc.  Each  experiment  is  now 
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The  SSS  sites  in  this  region  are  identified  during  the  appropriate  months.  High  nutrient  concentrations 
in  northern  latitudes  and  increased  sunlight  in  summer  months  result  in  an  intensive  high-latitude 
August  bloom. 


Figure  3  Mean  Monthly  Phytoplankton  Pigment  Concentration  in  the  North  Atlantic  Ocean  from  the 
Coastal  Zone  Color  Scanner 


-7- 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL,  MARYLAND 


CST-7 


CST-4 


February 


Phytoplankton 

Pigment 

Concentration 

April  (mg/m3) 


STB97-1 22-11 


Cloud  cover  limited  biological  productivity  during  CST-7  in  February  (upper  panei).  Spring  blooms 
likely  influenced  the  water  column  properties  during  CST-4  in  April  (lower  panel). 

Figure  4  Mean  Monthly  Phytoplankton  Pigment  Concentration  in  the  North  Pacific  Ocean  from  the 
Coastal  Zone  Color  Scanner 
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conveniently  described  by  the  set  of  mean  prediction  errors,  with  corresponding  base¬ 
line  environmental  descriptors,  listed  in  Table  2.  Note  that  wind  speed  values  do  not 
appear  in  Table  2.  The  objective  is  to  identify  additional  environmental  factors,  beyond 
local  wind  effects,  that  contribute  to  SSS  variability.  It  is  assumed  that  wind  speed  con¬ 
tributions  are  adequately  represented  by  the  ONE  model  and  are  hence  already 
'accoimted  for'  in  the  SSSe  values. 


Table  2 

Mean  Quantities  Used  for  Investigation  of  Site-to-Site  Differences 


Source 

h 

m 

sst 

°C 

chi 

O 

mg/m 

Mean  SSSe  (dB) 

0  - 150 

Hz 

400  -  600 
Hz 

800-1000 

Hz 

CST  1 

1.8 

10.7 

1.140 

1.1 

-1.3 

-  no  data  - 

CST2 

2.5 

8.3 

0.516 

4.8 

2.9 

3.2 

CST  3 

1.4 

28.1 

0.045 

-2.3 

-5.3 

-0.1 

CST  4 

3.2 

3.3 

0.320 

2.1 

-1.4 

0.8 

CST  5 

0.6 

23.0 

0.088 

2.6 

2.9 

5.0 

CST  7 

3.3 

5.5 

0.176 

1.0 

-0.7 

0.3 

Chapman- 

2.7 

18.0 

0.263 

-1.5 

-3.2 

-2.6 

Harris  * 

*  Chapman-Harris  hg  and  sst  estimates  obtained  from  U.S.  Navy  Marine  Climatic  Atlas  CD-ROM. 


The  combined  SSSe  values  depend  on  the  three  environmental  descriptors  in  a 
manner  that  agrees  with  physical  intuition.  To  facilitate  display  of  these  multidimen¬ 
sional  results,  we  will  first  normalize  the  data  with  respect  to  the  dependence.  Note 
that  CST-2  and  CST-4  are  both  cold  water  experiments  with  close  chi  values,  and  that 
CST-5  and  CST-3  are  both  warm  water  experiments  with  close  chi  values.  The  SSSe  dif¬ 
ferences  between  these  experiment  pairs  should  be  dominated  by  effects.  This  is  veri¬ 
fied  by  the  SSSe  vs.  plot  of  Figure  5.  The  differences  between  SSSe  values  of  each 
experiment  pair,  and  at  each  frequency  range,  are  represented  by  a  series  of  linear 
regressions.  Note  that  all  of  the  regressions  are  of  similar  slope  with  higher  observed 
scattering  strengths  (relative  to  ONE  predictions)  as  mean  wave  height  increases.  The 
regression  slopes  were  averaged  within  frequency  bands  and  used  to  produce  wave- 
height  normalized  prediction  errors  {SSSn)  for  each  of  the  observations  in  Table  2.  The 
normalization  is  given  by 


SSSn  =  m{h^i^gp-h^)  +  SSSe , 
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0-150  Hz  400-600  Hz  800-1000  Hz 


Mean  regression  slopes  for  each  frequency  band  were  used  to  normalize  the  SSSe  averages  in 
Table  2  for  surface  wave  effects. 


Figure  5  Influence  of  Background  Wave  Energy  on  Surface  Scatter  Strength  for  Sites  with  Similar 

Temperature  and  Chlorophyll  Characteristics 

where  m  is  the  average  regression  slope  for  each  frequency  band  indicated  on  Figure  5 
and  a  reference  wave  height  of  =  2.0  m  was  chosen. 

There  are  now  two  remaining  variables  on  which  SSSn  depends:  phytoplankton 
pigment  concentration  and  ocean  temperature.  Inspection  of  the  site  environmental  data 
in  Table  2  suggests  that  data  can  be  sorted  into  two  groups:  (1)  a  set  of  cold-water  obser¬ 
vations  (ssf  <  11  °C)  that  includes  CST-1,  CST-2,  CST-4,  and  CST-7;  and  (2)  a  set  of  low- 
productivity  observations  (chi  <  0.3  mg/m^)  that  includes  CST-3,  CST-5,  CST-7,  and 
Chapman-Harris.  Note  that  only  CST-7  falls  into  both  groups. 

The  dependence  of  SSSn  on  pigment  concentration  for  the  cold-water  observation 
set  appears  in  Figure  6.  The  results  imply  an  important  role  of  biological  activity  in  sur¬ 
face  scatter.  This  trend  is  strongest  at  mid-frequencies  (400-600  Hz).  Scattering  level 
increases  (represented  by  decreasing  SSSn)  with  biological  production  are  probably  due 
to  higher  dissolved  gas  levels  from  biological  productivity  and  the  presence  of  biologi¬ 
cal  surfactants,  both  of  which  will  extend  the  life  of  ambient  bubbles. 
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Influence  of  Biological  Productivity  on  /?3-Normalized  Surface  Scatter  Strength  Prediction 
Error  for  the  Cold  Water  {<  1 1  °C)  Observation  Sites 
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The  dependence  of  SSSn  on  sst  for  the  low  pigment  observations  appears  in  Fig¬ 
ure  7.  A  definitive  trend  of  increasing  scatter  (represented  by  decreasing  SSSn)  with 
increasing  temperature  is  observed  at  all  frequencies.  These  results  indicate  an  impor¬ 
tant  role  of  ambient  ocean  temperature  in  modulating  surface  scatter  levels.  The  influ¬ 
ence  of  temperature  on  bubble  entrainment  and  gas  dissolution  results  in  high  levels  of 
ambient  bubbles  in  warm  water.  Furthermore,  near-surface  stratification  in  warmer 
regions  may  effectively  trap  bubbles  near  the  surface  where  the  scattering  effect  is  most 
important. 

The  results  shown  in  Figures  5  through  7  suggest  that  a  linear  model  for  the 
dependence  of  SSSe  on  h^,  sst,  and  chi  might  explain  a  significant  fraction  of  the  site-to- 
site  variance  in  SSS.  The  method  of  least-squares  multiple  regression  was  employed  to 
test  the  performance  of  a  linear  model  for  the  surface  scatter  prediction  error 

SSSe  =  ah^  +  ^sst  +  ychl  +  b, 

with  fit  coefficients  a,  p,  and  y  for  wave  height,  ocean  temperature,  and  chlorophyll 
concentration,  respectively.  For  the  model  calculations,  the  SSSe  averaging  was 
restricted  to  narrow  wind  speed,  frequency,  and  grazing  angle  bands  so  that  the  results 
represent  a  specific  set  of  environmental  and  acoustic  conditions. 

This  preliminary  model  is  quite  successful  in  describing  gross  site-to-site  differ¬ 
ences  in  SSS.  Typical  results  appear  in  Figure  8.  Here  the  regression  coefficient  and  the 
three  model  fit  coefficients  are  plotted  as  a  function  of  acoustic  frequency  at  15-25°  graz¬ 
ing  angle  for  the  case  of  8-10  m/s  winds.  Note  that  a  high  regression  coefficient  is 
obtained  at  all  frequencies  indicating  that  the  model  explains  a  large  percentage  of  the 
total  variance.  There  is  a  linear  transition  of  all  model  parameters  across  frequency  with 
the  and  sst  fit  coefficients  essentially  constants. 

Using  a  and  P  to  normalize  the  data  for  both  wave  height  and  ocean  temperature 
effects,  as  demonstrated  earlier,  allows  a  direct  comparison  of  prediction  error  with  phy¬ 
toplankton  pigment  concentration.  This  normalization  is  given  by 

where  the  reference  values  =  2.0  m  and  =  15°C  were  chosen.  The  normal¬ 

ized  500  Hz  SSS  prediction  errors  (SSSn)  at  each  site  for  winds  of  ~10  m/s  appear  in 
Figure  9.  A  logarithmic  fit  to  the  data  is  given  by 

SSSn  =  -[4.8  +  3.61og(c/i/)], 
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Figure  8  Multiple  Linear  Regression  Results  for  the  Case  of  Uiq=  8-10  m/s  and  15-25°  Acoustic 
Grazing  Angle 

with  regression  coefficient  r  =  0.91.  The  CST-2  results  are  slightly  farther  than  two  stan¬ 
dard  deviations  from  the  mean  (see  Figure  2)  and  are  not  included  in  the  regression. 
This  is  perhaps  due  to  inaccuracy  of  the  monthly  mean  CZCS  chi  value  in  representing 
CST-2  conditions^  inspection  of  the  CZCS  map  of  Figure  3  (middle  panel)  suggests 
extreme  spatial  variability  in  chi  at  that  location  and  time.  The  remaining  data  of 
Figure  9  indicate  that  phytoplankton  blooms  contributed  up  to  5  dB  in  SSS  site-to-site 
differences  during  CST.  Extrapolating  these  results  to  nearshore  phytoplankton  pigment 
concentrations  of  50  to  100  mg/m^  indicates  that  SSS  values  can  be  10  to  12  dB  higher 
than  ONE  model  predictions  in  shallow  water  environments. 
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Phytoplankton  Pigment  (mg/m^) 

Extrapolation  to  nearshore  conditions  with  chi  ~  50-100  mg/m  indicates  ONE  model  SSS 
predictions  can  be  -  10-12  dB  too  low  in  highly  productive  shallow  water  environments. 

Figure  9  Dependence  of  Normalized  SSS  Prediction  Error  on  Phytoplankton  Pigment  Concentration 
for  an  Acoustic  Frequency  of  500  Hz  at  10  m/s  winds,  20°  Grazing  Angle  and  with  hg  =  2.0 
m  and  sst=  15°C 

5.  Surface  Wave  Descriptors  for  Bubbles 

The  CST-7  experiment,  conducted  during  February  1992  in  the  dynamic  Gulf  of 
Alaska,  provided  an  excellent  wind,  wave,  and  whitecap  observation  set  on  which  to 
test  the  hypothesis  that  surface  wave  parameters  can  be  used  to  replace  or  augment 
wind  speed  descriptors  in  bubble-related  air-sea  process  models  {Reference  3).  Wave 
spectral  parameters  were  obtained  from  a  newly  developed  spectral  partitioning  and 
swell  tracking  approach.  Using  these  parameters,  the  relationships  between  wind  forc¬ 
ing  variability,  wind  sea  growth,  swell,  and  wave  breaking  in  the  open  ocean  were 
investigated. 
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5.1  Automated  Wave  Spectral  Partitioning,  Swell  Tracking,  and  Storm 
Source  Identification 

To  facilitate  extraction  of  wave  descriptors  from  directional  wave  spectra,  we 
have  developed  a  fully  automated  method  for  wave  spectral  partitioning,  swell  track¬ 
ing,  and  storm  source  identification  using  MATLAB  programming  tools  (References  4  and 
5).  Following  an  approach  suggested  by  Susanne  Hasselmann  (Reference  6),  peaks  in 
wave  buoy  directional  spectra  are  isolated  based  on  topographic  minima,  with  wind  sea 
peaks  identified  by  wave  age  criteria.  We  have  added  a  swell  tracking  algorithm  that 
provides  a  unique  approach  to  storm  source  prediction  when  combined  with  linear 
wave  theory.  Until  now,  identification  of  swell  source  times  and  locations  from  wave 
spectra  was  a  laborious  manual  procedure  and  limited  to  the  dominant  swell  system. 

Over  the  6-day  CST-7  wave  record,  44  swell  systems  were  identified  with  up  to 
three  systems  coexisting  at  any  given  time  (Reference  7).  The  presence  of  atmospheric 
disturbances  on  surface  weather  charts  validated  the  storm  source  predictions  for  over 
85  /o  of  these  systems.  A  synthesis  of  results  from  the  swell  and  wind  sea  separations, 
appearing  in  Figure  10,  provides  a  distinctive  representation  of  the  surface  wave  history 
over  the  duration  of  the  experiment. 

The  wave  partitioning  technique  has  several  important  applications.  It  would 
facilitate  the  reduction  of  large  wave  data  sets  such  as  those  from  satellite  observations 
and  global  wave  evolution  models.  Furthermore,  the  automated  swell  tracking  and 
source  identification  procedure  has  a  potentially  important  application  in  the  forecast¬ 
ing  of  storm  waves  for  shipping  routes  and  populated  coastal  regions.  In  the  following 
section,  we  describe  how  the  results  are  used  to  evaluate  the  use  of  surface  wave 
descriptors  for  process  models  related  to  near-surface  bubble  generation. 

5.2  Surface  Wave  Dissipation  and  Whitecap  Formation 

Wave  breaking  is  involved  in  many  important  air-sea  processes,  including  wave 
energy  dissipation;  momentum,  gas,  and  heat  fluxes;  ambient  noise  generation  by  air 
entrainment;  radar  backscatter  by  whitecaps;  and  acoustic  reverberation  from  bubble 
clouds.  However,  the  mechanisms  of  wave  breaking  in  the  ocean  are  not  well  under¬ 
stood;  as  a  result,  many  of  these  related  processes  are  empirically  modeled  using  wind 
speed  parameters  (Reference  8).  We  have  found  that  wave  spectral  parameters  can  be 
used  in  place  of  wind  parameters  in  the  estimation  of  processes  related  to  wave 
breaking. 
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Lower  panel  wind  vectors  indicate  the  direction  towards  which  the  air  is  moving.  The  upper 
panel  depicts  the  evolution  of  the  principal  wind  sea  and  swell  wave  systems  over  the  obser¬ 
vation  period.  Wave  height,  direction,  and  frequency  are  represented  by  vector  length,  azi¬ 
muth,  and  ordinate  position,  respectively. 


Figure  10  Wind  and  Wave  Vector  History  from  CST-7  Gulf  of  Alaska  Observations 


The  total  rate  of  wave  energy  dissipation,  e, ,  was  estimated  from  the  CST-7  wind 
sea  spectral  partitions  using  equilibrium  range  concepts  {Reference  9).  As  demonstrated 
in  Figure  11,  replacing  wind  speed  with  dissipation  rate  in  the  standard  power-law 
description  of  whitecap  fraction  W  resulted  in  a  new  proportionality, 

1.5 

W  °<=  Ej  , 

which  reduces  the  scatter  of  W  about  the  independent  variable  by  2  to  3  orders  of  mag¬ 
nitude.  These  improvements  result  primarily  from  the  removal  of  wind  trend  effects 
that  add  uncertainty  to  wind  speed  relationships.  The  successful  description  of  white- 
caps  by  our  dissipation  rate  estimates  lends  credibility  to  the  equilibrium  range  theory 
as  a  viable  approach  to  modeling  bubble  generation  by  wave  breaking. 
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Wave  dissipation  rate  reduces  the  scatter  about  a  power-law  model. 

Figure  1 1  Dependence  of  Whitecap  Fraction  on  Wind  Speed  and  Wave  Dissipation  Rate 

Remote  sensing  techniques,  wave  buoy  arrays,  and  regional  wave  modeling 
efforts  are  continually  being  refined  so  that  soon  high-quality  directional  wave  spectra 
will  be  routinely  available  on  a  global  scale.  It  is  conceivable  that  air-sea  process  models 
for  such  mechanisms  as  gas  flux,  radar  backscatter,  underwater  ambient  noise,  and  near¬ 
surface  acoustic  reverberation  will  soon  improve  significantly  due  to  the  use  of  surface 
wave  spectral  parameters  in  place  of  external  forcing  parameters.  Since  wave-related 
processes  are  dynamically  linked  to  intrinsic  properties  of  the  wave  field,  the  use  of 
wave  parameters  should  improve  model  performance. 

6.  Conclusions  and  Recommendations 
6.1  Conclusions 

•  Gross  site-to-site  differences  in  SSS  are  explained  by  seasonal  and  geo¬ 
graphic  environmental  factors. 

•  Biological  productivity  significantly  increases  SSS  in  the  open  ocean.  This 
is  likely  due  to  increased  dissolved  gas  levels  and  the  presence  of  biochem¬ 
ical  surfactants,  both  of  which  will  extend  bubble  Ufetimes.  This  effect  will 
be  amplified  in  shallow  water  environments. 
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•  SSS  is  higher  in  warm  water.  This  is  likely  a  result  of  greater  bubble  supply 
due  to  viscosity  and  surface  tension  effects  and  decreased  gas  solubility 
which  acts  to  extend  bubble  lifetimes.  Increased  temperature  variability 
will  greatly  influence  scattering  strength  statistics  in  shallow  water  envi¬ 
ronments. 

•  SSS  increases  in  higher-energy  environments.  Here  the  background  ocean 
energy  level  was  characterized  by  significant  wave  height  averaged  over 
several  days. 

•  Wave  spectral  partitioning  can  be  employed  to  extract  wave  descriptors  for 
bubble-related  processes  in  the  ocean. 

•  Wave  dissipation  rate  improves  the  standard  power-law  model  for  the  pre¬ 
diction  of  ocean  whitecaps. 

6.2  Recommendations  for  Continued  Research 

•  Prepare  a  robust  SSS  model  with  physical  and  biological  parameters  for 
site-to-site  variability 

•  Determine  the  relevance  to  other  SSS  observation  sets  including  ASREX 

•  Estimate  SSS  statistics  for  various  shallow  water  strategic  areas  of  interest 
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